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^ ; ABSTRACT 

fC ■ We present new results on the Lya emission-line kinematics ofl8z~2 — 3 star- forming galaxies 

?pH| with multiple-peaked Lya profiles. With our large spectroscopic database of UV-selected star-forming 

Q , galaxies at these redshifts, we have determined that ~30% of such objects with detectable Lya emis- 

^ ■ sion display multiple-peaked emission profiles. These profiles provide additional constraints on the 

^ ■ escape of Lya photons due to the rich velocity structure in the emergent line. Despite recent advances 

Ci . in modeling the escape of Lya from star- forming galaxies at high redshifts, comparisons between mod- 

els and data are often missing crucial observational information. Using Keck II NIRSPEC spectra 
of Ra {z ~ 2) and [OIIIJA5007 {z ^ 3), we have measured accurate systemic redshifts, rest-frame 
optical nebular velocity dispersions and emission-line fluxes for the objects in the sample. In addition, 
rest-frame UV luminosities and colors provide estimates of star-formation rates (SFRs) and the degree 
of dust extinction. In concert with the profile sub-structure, these measurements provide critical con- 
straints on the geometry and kinematics of interstellar gas in high-redshift galaxies. Accurate systemic 
^-u . redshifts allow us to translate the multiple-peaked Lya profiles into velocity space, revealing that the 

• ■ majority (11/18) display double-peaked emission straddling the velocity- field zeropoint with stronger 

^^ \ red-side emission. Interstellar absorption-line kinematics suggest the presence of large-scale outflows 

^^ . for the majority of objects in our sample, with an average measured interstellar absorption velocity 

offset of (Avabs) = —230 km s~^. A comparison of the interstellar absorption kinematics for objects 
with multiple- and single-peaked Lya profiles indicate that the multiple-peaked objects are charac- 
terized by significantly narrower absorption line widths. We compare our data with the predictions 
of simple models for outflowing and infalling gas distributions around high-redshift galaxies. While 
popular "shell" models provide a qualitative match with many of the observations of Lya emission, we 
find that in detail there are important discrepancies between the models and data, as well as problems 
with applying the framework of an expanding thin shell of gas to explain high-redshift galaxy spectra. 
Our data highlight these inconsistencies, as well as illuminating critical elements for success in future 
models of outflow and infall in high-redshift galaxies. 

Subject headings: galaxies: high-redshift — galaxies: formation — galaxies: ISM — radiative transfer 
— line: profiles 

1. INTRODUCTION Feedback commonly refers to large-scale outflows of 

The process described as "feedback" is considered "^fss metals, energy, and momentum from galaxies, 

a crucial component in models of galaxy formation. "^^'f" therefore regulate the amount of gas available 

to lorm stars, as well as the thermodynamics and 

chemical enrichment of the surrounding intergalactic 

^ Based, in part, on data obtained at the W.M. Keck Obser- medium (IGM). The evidence for feedback in high- 

yatory, which is operated as a scientific partnership a^^^ redshift star-forming galaxies comes in several forms, 
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l2QlQf ) , and the nature of the IGM environments of vigor- 
ously star- forming galaxies, in terms of the optical depth 
and kinemati cs of the surrounding H I and heavy ele- 
ments (Adelb erger et ah 2003, 2005). Yet, in spite of ev- 
idence for ubiquitous outflows at high redshift, estimates 
of fundamental physical quantities such as gas column 
densities and mass outflow rates in the superwinds have 
remained elusive. While there is much observational evi- 
dence for outflows in high-redshift galaxies, both analytic 
models and hydrodynamical simulations of these systems 
suggest that, at the same cosmic epochs, they should be 
rapidly accreting cold gas fro m the IGM, fueling their 
active rates of star formation (jBirnboini fc Dekell 120031 : 
iKeres et al.ll2005L 120091 : iDekel et al.ll2009f ). Searching for 
observational signatures of the process of cold gas accre- 
tion at high redshift remains an open challenge. 

The Lya feature is one of the most widely used probes 
of star formation in both the nearby and very distant 
universe. While Lya photons are initially produced by 
recombining ionized gas in H II regions, resonant scat- 
tering through the interstellar medium (ISM) of galaxies 
can lead to extreme modulation of the intrinsic emis- 
sion profile, in both frequency and spatial location. Ab- 
sorption by dust can completely suppress Lya emission, 
producing a strong absorption profile, even in galaxies 
with high rates of star formation. It is therefore dif- 
ficult to determine the intrinsic properties of the gas 
giving rise to Lya emission from the observed profile 
alone. Only recently has it become possible to make de- 
tailed theoretical predictions for the Lya profiles emer- 
gent from con iplex systems similar to those observed at 
hi^h redshift (lAhn et al."2002': "Zh eng fc Mir alda-Escude* 
120021 : iHansen fc O h 2006; Verham me et al . 2006, 2008). 
These new calculations (e.g., Verhamme et al. 2006) em- 
ploy a Monte Carlo approach to propagate a representa- 
tive ensemble of Lya photons through gas and dust of ar- 
bitrary spatial and velocity distribution outputting Lya 
profiles as would be observed. By comparison with Lya 
profiles in actual galaxy spectra, it is possible, in prin- 
ciple, to recover quantities such as the expansion/infall 
velocity of the outflow/inflow, the column density and ve- 
locity dispersion of absorbing gas, and the gas covering 
fraction. Therefore, modeling of observed Lya emission 
profiles represents an independent method of probing the 
processes of feedback and accretion at high redshift. 

The majority of Lya emission profiles at high red- 
shift fall in the category of single-peaked and asymmetric 
(jShaplev et al .112 0031 : Ifapken et al.ll2007p. which is a nat- 
ural outcome of an expanding medium () Verhamme et al] 
I2OO8). However, multiple-peaked Lya profiles, seen in a 
fraction of star- forming UV-selected galaxies at z ~ 2 — 3, 
offer a particularly detailed perspective on Lya escape 
due to the rich structure of the emergent line. In prin- 
ciple, the structure of the Lya line encodes the velocity 
field and density distribution of the gas through which it 
has emerged. For example, a symmetric double-peaked 
profile centered on the velocity-field zeropoint is a natural 
outcome of the radiative transfer of Lya photons through 
a static medium ( Osterbrock 1962 ). Recent advances 
in modeling the escape of the Lya photons from star- 
forming galaxies at 2; ~ 2 — 3 have isolated several fea- 
tures of the Lya profile that may be ex pected for specific 
gas geometries and velocity fields (e.g. 'Verha mme et al.l 
12006^ JLaursen et ah , 2QD9%b; Barnes et a l, 20111) . Since 



Lya is the most readily observed high-redshift emission 
line, there is enormous potential to better understand 
the structure of the early galaxies by interpreting Lya 
line profiles. 

Previous attempts to compare observed Lya line mor- 
phologies to theoretical predictions have been missing 
crucial information, weakening the derived constraints 
on galaxy outflow and inflow properties. In particular, 
accurate systemic redshift measurements, nebular line 
widths, and intrinsic ionizing photon fluxes have been 
absent fro m mo st previous comparisons (but see, e.g., 
lYang et al.lf2011l : lMcLinden et al.ll2011f ). These three ob- 
servables are critical for anchoring the velocity scale of 
the models, constraining the mass and thermal motions 
of the gas, and determining the overall normalization for 
a given model. In this paper, we use new results ob- 
tained from Ha and [OIIIJA5007 emission lines to provide 
the critical missing constraints on the observed kinemat- 
ics of star-forming galaxies at 2: ~ 2 — 3 with multiple- 
peaked Lya emission. In addition, we have used our large 
database of Lya emission lines in high-redshift objects to 
determine the frequency of the multiple-peaked systems, 
providing a global perspective on the potential of using 
Lya morphology to reveal z ^ 2 — 3 gaseous structure. 

Our method of target selection from the parent sam- 
ple of UV-selected galaxies at z ~ 2 — 3 is explained 
in Section [21 while the observations and data reduction 
are described in Section [3l The Lya velocity profiles 
are presented in Section |4] with precise velocity- field ze- 
ropoints determined from our Ha (or [OIIIJA5007) mea- 
surements. In Section [5] the measured physical quanti- 
ties for each system are reported. Section [6] describes 
current Lya radiative transfer models, in addition to a 
qualitative comparison between some simple models and 
our measured Lya velocity profiles. Finally, in Section El 
we summarize our results and discuss how our measure- 
ments will be used in future work to accurately model 
the processes that give rise to the multiple-peaked Lya 
profiles from star- forming galaxies at z ~ 2 — 3. We as- 
sume a flat ACDM cosmology with Q^ = 0.3, Qa = 0.7, 
and Ho = 70 km s~^ Mpc~^ 

2. SAMPLE SELECTION 

Our target galaxie s were drawn f rom the UV -selected 
sample described in iSteidel et all (J2003L 120041 ). These 
UV-selected surveys have yielded over 3,000 spectroscop- 
ically confirmed star- forming galaxies at z ~ 1.5 — 3.5. 
The majority of the spectra cover the HI Lya feature at 
1216 A. The observed Lya profiles for UV-selected star- 
forming galaxies vary widely in strength, from damped 
absorption to strong emission. Of the ~1500 objects that 
show net emission, a fraction exhibit a multiple-peaked 
profile. While this phenomenon has been previ ously re- 
ported (jTapken et al.l 120071 : iQuider et all l2QQ9f ) its fre- 
quency of occurrence has never been systematically ana- 
lyzed. We have used our database of star-forming galax- 
ies to assess the frequency of multiple-peaked Lya pro- 
files in z '^ 2 — 3 galaxies in order to determine whether 
such profiles represent rare outliers or are commonplace 
features of the galaxy population. 

To identify the frequency of multiple-peaked profiles, 
we first separated the full spectroscopic sample into 
subsets of galaxies with and without detectable Lya 
emission. The ~1500 galaxies with Lya emission were 
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Fig. 1. — 7^-band images (rest-frame UV) for our sample of 18 objects. In each image, the NIRSPEC long slit position is overlaid with 
blue, solid lines. The LRIS slit is overlaid with red, dashed lines. All images are 12^^x12^^, with north up and east to the left. Typical 
seeing during these exposures was O.^'^S — 1/^0. 



then considered for further study. Ah of the optical 
(rest-frame UV) spectra were obtaine d usin^ the LRI S 
spectrograph on the Keck I telescope (jOke et al.lll995f ). 
These spectra were taken with the 300 line mm~^ grat- 
i ng blazed at 5000 A, or, following the LRIS-B upgrade 
(jSteidel et alJl2004D . with the 400 line mm-^ or 600 line 
mm~^ grism blazed at 3400 A, and 4000 A, respectively. 
The resolution of the spectra (taken through l.'"2 slits) 
establishes a lower limit in velocity space on the separa- 
tion of multiple peaks that can be identified. From this 
sample we were able to identify objects that have a min- 
imum peak separation of 225 [370, 500] km s~^ for the 
600 [400,300]-line grism/grating. We separated the spec- 
tra by grism/grating and analyzed each resolution group 
individually. Approximately 60% of the spectra were ob- 
tained with the 400- and 600-line grisms, while the rest 
of the spectra were taken using the 300-line grating. 

We used a systematic algorithm to search for multiple 
peaks within the Lya profile. Several criteria were used 
to minimize contamination from noisy spectra which can 
yield a false indication of multiple peaks due to noise 
spikes. The first criterion was a signal-to-noise (S/N) 
lower limit of 3 on the strongest fiux peak height, which 
was implemented to establish a consistent threshold for 
characterizing Lya profiles. This criterion reduced our 
sample to ~1000 spectra with 163, 527, and 402 spec- 
tra, respectively, from the 600- and 400-line grism, and 
300-line grating. The spectra were then separated into 
resolution groups and the algorithm was implemented 
independently on each sub-sample. First, we searched 
between 1210 and 1225 A in each spectrum for the max- 
imum fiux value, given that this wavelength range typ- 
ically contains the observed Lya peak in the spectra of 
star- forming galaxies. This maximum was identified as 
the primary Lya emission peak. The continuum level 



was then estimated for the blue and red side separately 
from the mean of the fiux over a ~100 A span blueward 
and redward of 1210 and 1225 A, respectively. Starting 
from the primary peak's maximum fiux point, we stepped 
down towards the blue side until reaching the continuum 
level. If a local minimum was noted before reaching the 
continuum we marked the local maximum after that min- 
imum as another possible peak for the Lya line. Another 
criterion used to avoid noise spikes was for the secondary 
maximum to span at least two adjacent steps in wave- 
length. This procedure was then repeated on the red side 
of the primary peak. If an additional peak was noted by 
the algorithm, we required the fiux ratio of the additional 
peak to the primary peak to be greater than -^ , with N 
the S/N of the primary peak. This value was chosen in 
order to reject any small noise spikes in the identification 
of the additional peak. Our sample is therefore complete 
for primary to secondary ratios less than or equal to 3. At 
higher primary-peak S/N ratio we can also probe larger 
ratios of primary-to-secondary peak heights. In practice, 
however, the typical primary-to-secondary peak ratio is 
~ 3. In addition to analyzing the full sample of spec- 
tra with our objective criteria, we examined the spectra 
by eye for any objects that might have been missed by 
the algorithm described above. The additional multiple- 
peaked objects found by eye comprise a small percentage 
of the total sample, roughly ^ 15%. 

For the 300-line grating, the percentage of spectra with 
multiple-peaked Lya emission identified is 23% (92 spec- 
tra). For the 600- and 400-line grisms the values are, 
respectively, 27% (44 spectra) and 20% (103 spectra). 
In addition to the main spectroscopic sample, a set of 
121 objects was followed up with 400-line grism spectra 
using significantly longer exposure times (7 — 13 hours 
as opposed to 1.5 hours), with correspondingly higher 
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S/N (Bogosavljevic 2010, Ph.D Thesis). In this sample, 
97 objects show Lya emission, and 33% (32 spectra) of 
the emitters are identified as containing multiple-peaked 
Lya emission using the criteria described above. Higher 
resolution (i.e., 600-line grism) and S/N (i.e., deeper 400- 
line grism) spectra are more likely to have identifiable 
multiple-peaked Lya emission if it is present. Based on 
the 600-line and deep 400-line samples we therefore assert 
that the prevalence of multiple-peaked profiles among 
objects with Lya emission is ~30% in a sample of UV- 
selected star- forming galaxies at z ~ 2 — 3. This percent- 
age references a S/N lower limit of 3 on the strongest fiux 
peak and a resolution limit of 225 km s~^. 

In order to access the rest-frame optical nebular emis- 
sion lines (Ra at 2; ^ 2 and [OIIIJA5007 at z ~ 3) en- 
abling the determination of systemic velocities for our 
galaxies, we are restricted to those objects in the full 
multiple-peaked sample whose redshifts ensure that ei- 
ther Ha or [OIIIJA5007 falls into a window of atmo- 
spheric transmission. The sharply increasing thermal 
background past 2.3 jam further limits the redshift range 
over which rest-frame optical emission lines are acces- 
sible. In the iT-band the available redshift range is 
z - 2.0 - 2.6 for Ra and z - 2.9 - 3.7 for [OIII]. In the 
iJ-band [OIII] can be accessed from z ~ 1.9 — 2.6. After 
applying the restricted redshift ranges to the multiple- 
peaked objects determined from the algorithm described 
above, we identified 193 objects as potential targets for 
near-infrared spectroscopic follow-up. From this sample 
we selected the 18 objects described in this paper, which 
are representative in terms of their Lya kinematics (see 
Section HI . 

3. OBSERVATIONS AND DATA REDUCTION 

Both near-IR and optical spectra are required for our 
analysis of the kinematics of Lya emission profiles. In 
this section, we describe both types of spectroscopic data, 
as well as our basic empirical measurements and uncer- 
tainties. 

3.1. NIRSPEC Observations and Data Reduction 

Near-IR spectra were obtain ed using the NIRSPEC in- 
strument ([McLean et al.|[T998[ ) on the Keck II telescope. 
Our data were collected during three observing runs in 
2009 July, August, and October, for a total of 5 nights. 
On average the total exposure time for each object was 
4x900 seconds, though the number of 900 second ex- 
posures per target ranged from 3 to 6 for a total inte- 
gration time of 0.75 to 1.5 hours. All targets were ob- 
served with a 0.''76 x 42'' long slit. For most objects 
we used the N6 filter, a broad H -\- K filter centered at 
1.925/im with a bandwidth of 0.75/im. For objects at 
2;^ 2.5 — 2.6, we used the N7 filter, which is centered at 
2.23/im with a bandwidth of 0.80/im. The spectral res- 
olution as determined from sky lines was ~ 15 A, which 
corresponds to Av -230 km s"^ (R:^1300) for the N6 
filter and A^; -200 km s"^ (R::^1500) for the N7 filter. 
Conditions were photometric, with the seeing ranging 
from 0.''33 - 0.'70. 

Due to the faint nature of these objects in the iT-band 
we acquired each target using blind offsets from a bright 
star in the surrounding field. We returned to the off- 
set star between each integration of the science target 



to recenter and dither along the slit. In nearly all cases 
we dithered back and forth between two positions near 
the center of the slit. The position angle was chosen to 
match that of the LRIS observation unless the object 
appeared extended at a different angle in the existing 
7^-band images (S teidel et al . 2003, 2004). We measured 
Ha or [OIIIJA5007 for a total of 18 objects during our 
observing runs, which were drawn from the parent sam- 
ple discussed in Section [2l All observational information 
can be found in Table [TJ 7^-band images of the objects 
with the NIRSPEC and LRIS slits overlaid are presented 
in Figure [H 

Data red uction was performed following the method 
described in lLiu etldl |2008), where the sky background 
was subtracted from the two-dimensional unrectified sci- 
ence images using an optimal method (Kelson 2003, G. 
D. Becker 2006, private communication). After back- 
ground subtraction, cosmic rays and bad pixels were re- 
moved from each exposure, which was then rotated, cut 
out along the slit, and rectified in two dimensions to 
take out the curvature both in the wavelength and spa- 
tial directions. The sky lines were fit with a low-order 
polynomial and a b-spline fit was used in the disper- 
sion direction. The final rectified two-dimensional expo- 
sures for each object were then registered and combined 
into one spectrum. A one-dimensional spectrum was ex- 
tracted from the two-dimensional reduced image along 
with the corresponding la error spectrum. The aver- 
age aperture size along the slit was 2.'"1, with a range of 
l.'"7 to 2. ''9. The spectrum was then flux-calibrated us- 
ing A-type stars observed during o ur NIRSPEC runs ac- 
cord ing to the metho d described in lShaplev et al.l (|2005f ) 
and lErb et al.l (|2003f ). Finally the flux- calibrated, one- 
dimensional spectrum was placed in a vacuum, heliocen- 
tric frame. 

3.1.1. Objects With Spatially-Resolved NIRSPEC Spectra 

While the majority of the objects in our sample have 
spatially unresolved Ha (or [OIIIJA5007) emission-line 
spectra, three objects exhibit spatially-extended nebu- 
lar emission: Q0449-BX167, Q2206-BX151, and Q1549- 
M22. If spatially-extended emission is coupled with rota- 
tion or velocity shear, the placement or position angle of 
the long slit has the potential to yield nebular emission 
centroid and velocity dispersion values that differ from 
the true properties of the galaxy. The differences may re- 
sult if the slit does not evenly sample both the approach- 
ing and receding parts of the rotation curve or if sys- 
tematic shear across the galaxy is interpreted as purely 
random motions. The three objects listed above were 
analyzed in more detail to determine whether their com- 
plexity might have affected our measurements. In the 
7^-band image for Q0449-BX167 there are two distinct 
components and our slit PA was chosen to have both 
components included. Q0449-BX167 has two separate 
components in the two-dimensional NIRSPEC spectrum 
that are offset in the spatial direction by 1.''3, but not 
in the spectral direction. In the 7^-band image, Q2206- 
BX151 is extended and our slit PA was chosen to encom- 
pass the extended region. In the two-dimensional NIR- 
SPEC spectrum Q2206-BX151 is extended in the spa- 
tial direction for a total of 1.'''3. There is a small tail 
that extends off of the main nebular emission, which is 
slightly offset in the spectral direction. The flux from 
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TABLE 1 
Observation Log with Keck II NIRSPEC 



Object 


R.A. (J2000) 


Dec. (J2000) 


^HIl"'' 


7^(mag) 


Exposure (s) 


Filter 


Date (UT) 


Q0100-MD40 


1 03 18.971 


13 17 02.841 


2.2506 


24.35 


4 


X 


900 


N6 


Aug. 


09, 2009 


Q0100-BX288 


1 03 20.931 


13 16 23.569 


2.1032 


23.75 


2 


X 


900 


N6 


Aug. 


09, 2009 


Q0449-BX167 


4 52 24.624 


-16 39 15.284 


2.3557 


24.02 


4 


X 


900 


N6 


Oct. 


10, 2009 


Westphal-DD29 


14 17 25.936 


52 29 32.219 


3.2401 


24.82 


4 


X 


900 


N6 


Aug. 


10, 2009 


Westphal-BX154 


14 17 32.173 


52 25 51.044 


2.5954 


23.96 


3 


X 


900 


N7 


July 


02, 2009 


Q1549-BX102 


15 51 55.977 


19 12 44.220 


2.1932 


24.36 


3 


X 


900 


N6 


July 


01, 2009 


Q1549-C20 


15 52 00.396 


19 08 40.773 


3.1174 


24.87 


3 


X 


900 


N6 


Aug. 


10, 2009 


Q1549-M22^ 


15 52 02.703 


19 09 40.011 


3.1535 


24.85 


1 


X 


900 


N6 


Jul. 


03, 2009 












2 


X 


900 


N6 


Aug. 


10, 2009 


Q1549-BX202 


15 52 05.090 


19 12 49.818 


2.4831 


24.37 


4 


X 


900 


N6 


July 


01, 2009 


Q1623-BX129 


16 25 28.732 


26 49 19.182 


2.3125 


24.04 


5 


X 


900 


N6 


July 


02, 2009 


Q2206-BX151 


22 08 48.650 


-19 42 25.569 


2.1974 


24.03 


4 


X 


900 


N6 


July 


01, 2009 


SSA22a-D17^ 


22 17 18.845 


18 16.667 


3.0851 


24.27 


6 


X 


900 


N6 


July 


02, 2009 


SSA22a-C49^ 


22 17 19.810 


18 18.729 


3.1531 


23.85 


6 


X 


900 


N6 


July 


02, 2009 


SSA22a-C31 


22 17 22.885 


16 09.492 


3.0176 


24.61 


3 


X 


900 


N6 


July 


03, 2009 


SSA22a-C32 


22 17 25.629 


16 13.095 


3.2926 


23.68 


4 


X 


900 


N6 


July 


03, 2009 


Q2343-BX70 


23 45 55.934 


12 44 50.317 


2.4086 


24.94 


4 


X 


900 


N6 


Aug. 


09, 2009 


Q2343-BX231 


23 46 20.108 


12 46 16.812 


2.4999 


25.15 


3 


X 


900 


N6 


Aug. 


10, 2009 


Q2346-BX181 


23 48 31.827 


21 39.077 


2.5429 


23.28 


3 


X 


900 


N7 


July 


01, 2009 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, 
arcminutes and arcseconds. 

^ Rest-frame optical nebular emission-line redshifts. Objects with z ~ 2.0 — 2.6 have Ho; measurements. Objects 
with z > 3 have [OIII]A5007 measurements. 
^ Typical uncertainty in zhii is A2; ~ 10~^. 

^ Object Q1549-M22 was observed on both July 3^ and August 10 and the science images from these nights 

were stacked together. 

^ Objects SSA22a-D17 and SSA22a-C49 were obtained on the same slit. 
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Fig. 2. — From left to right, the two-dimensional LRIS and NIR- 
SPEC spectra for Q0449-BX167, Q2206-BX151, and Q1549-M22. 
These three objects exhibit spatially-extended nebular emission. 
The horizontal and vertical axes represent, respectively, velocity 
with respect to the velocity zeropoint and spatial distance along 
the slit. In each image, the dashed vertical line represents the 
systemic velocity of the system {y = 0), and the horizontal axis 
spans ±1,500 km s"-*^ in velocity. The vertical axes span ±3.5^^ 
in the spatial direction, and a vertical scale bar of V is shown 
in the upper left-hand panels for Q0449-BX167. (Top) LRIS 2D 
Lyo; spectra. (Bottom) NIRSPEC 2D spectra. The spectra of 
Q0449-BX167 and Q2206-BX151 show the Ho; feature, while the 
spectrum of Q1549-M22 is centered on [OIII]A5007. Q0449-BX167 
has two components separated by l.''^3 along the slit, Q2206-BX151 
is extended spatially by l.''^3, and Q1549-M22 consists of two com- 
ponents separated by l.''^l in the spatial direction and 130 km s"-*^ 
in the spectral direction. 



the tail is comparable to the noise seen in the image 
and was not extracted. Since the structure in the two- 



dimensional spectra for Q0449-BX167 and Q2206-BX151 
(excluding the tail) is not shifted in the spectral direction 
there is no significant evidence for rotation or shear, and 
the nebular emission centroid and dispersion values for 
these objects should be robust. In fact, Q0449-BX167 
and Q2206-BX151 have among the smallest velocity dis- 
persion values of the sample. The third object, Q1549- 
M22, is slightly extended in the 7^-band image and again 
we chose our slit PA to correspond to the extended re- 
gion. There are two distinct components in the two- 
dimensional spectrum for Q1549-M22 and they are offset 
by l.''^! and 9 A, respectively, in the spatial and spectral 
directions. The two-dimensional NIRSPEC spectra for 
these three objects are presented, along with their corre- 
sponding two-dimensional LRIS spectra, in Figure [2l 

Object Q0449-BX167 was originally extracted to in- 
clude only the brighter component on the left hand side 
of the two-dimensional spectrum. To determine if the 
spatially-extended nebular emission affects our measure- 
ments, we re-extracted Q0449-BX167 to include both 
components seen in the two-dimensional spectrum. For 
Q0449-BX167 the centroid did not change and the veloc- 
ity dispersion changed by Acr^ = +2.5 km s~^ when com- 
paring the more extended extraction to our original one. 
This difference is within our errors. For our analysis we 
utilize the values obtained from our original extraction of 
Q0449-BX167. Object Q2206-BX151 was originally ex- 
tracted to include only the main elongated component. 
We re-extracted to include the tail for comparison. The 
centroid of the re-extraction did not change compared 
to the centroid measured from the original extraction. 
The velocity dispersion changed from an upper limit of 
41 km s~^ (not including the tail) to a measurement of 
63 km s~^ (including the tail). The re-extractions show 
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that the centroid measurements used in our analysis for 
Q0449-BX167 and Q2206-BX151 are robust. We include 
the velocity dispersion values measured from the origi- 
nal extraction for Q2206-BX151 in our analysis with the 
caveat that the added complexity might affect our mod- 
eling comparison (see Section [6]). 

Our fiducial measurements for Q1549-M22 are based 
on the extraction of both spectral components. Further 
analysis was conducted to understand how the spectral 
offset between the two components might affect our mea- 
surements. Each component was extracted and reduced 
separately for comparison with the integrated extraction 
of both components. Relative to the velocity zeropoint 
measured from the summed components, the offsets mea- 
sured for each component separately are A^' = 110 km 
s~^ and Av = —76 km s~^. The measured velocity 
dispersion of the spectrum integrated over both compo- 
nents (153 km s~^) is approximately equal to the sum 
in quadrature of the velocity dispersion of the individual 
components (125 km s~^ and 95 km s"-*^). The shift from 
the velocity zeropoint of each component and the large 
measured velocity dispersion from the summed spectra 
(the largest in our sample) may be evidence of two sepa- 
rate components giving rise to the two Lya peaks for this 
object. We decided to include Q1549-M22 in our sample 
with the caveat that its complexity in both [OIIIJA5007 
and Lya may signify a different underlying phenomenon 
giving rise to the multiple-peaked Lya emission from the 
one we will discuss in subsequent sections. 

3.2. LRIS Data Reduction 

The observations and reduction of the LRIS spectra 
for o ur targ:et sample ha ve been described in previous 
work CSteidel et al.l l2QQ3L[2QQl . These spectra were col- 
lected over the course of many years without consistent 
attention to the accuracy of the zeropoint of the wave- 
length solution at the 50 to 100 km s~^ level, given that 
their primary purpose was for the measurement of red- 
shifts and basic galaxy properties. In contrast, for the 
work presented here it is crucial to obtain a uniform and 
highly accurate wavelength zeropoint in order to place 
our LRIS and NIRSPEC data sets on the same veloc- 
ity scale. Accordingly, we re-calculated the wavelength 
solution of each LRIS spectrum. After computing the 
wavelength solution using the IRAF routine identify we 
adjusted the zeropoint of the solution such that bright, 
unblended sky lines (e.g., A5577.339, A6300.304) showed 
up at the correct wavelengths. The spectra were then 
corrected from an observed, air frame to a heliocentric, 
vacuum one. 

3.3. Line Centroid, Flux, & FWHM Measurements 

Measuring an accurate nebular emission centroid in or- 
der to obtain a precise systemic redshift for our objects 
was the principal objective of the NIRSPEC observations 
described in Section [3711 The flux and FWHM were also 
important measurements for our analysis. The lines that 
we set out to measure were Ha at 2; ~ 2 and [OIIIJA5007 
at z ^ 3. The centroid, flux, and FWHM were deter- 
mined by fitting a Gaussian profile to each emission line 
using the IRAF task, splot. 

A Monte Carlo approach was used to measure the un- 
certainties in the centroid, flux, and FWHM. For each ob- 
ject 500 fake spectra were created by perturbing the flux 



TABLE 2 
NIRSPEC Measurement Results 



Object 


2^HIl" 


FWHMhii^'^ 


i^HIl"'^ 


Q0100-MD40^ 


2.2506 


14.3 ± 1.0 


3.6 ± 0.3 


Q0100-BX288 


2.1032 


19.6 ± 0.4 


6.7 ± 0.1 


Q0449-BX167 


2.3557 


17.3 ± 1.4 


7.8 ± 0.7 


Westphal-DD29 


3.2401 


22.9 ± 1.3 


22.9 ± 0.5 


Westphal-BX154 


2.5954 


17.7 ± 2.0 


7.7 ± 0.8 


Q1549-BX102 


2.1932 


18.2 ± 1.1 


6.1 ± 0.3 


Q1549-C20 


3.1174 


17.3 ± 1.0 


7.3 ± 0.4 


Q1549-M22 


3.1535 


29.2 ± 2.8 


8.1 ± 0.6 


Q1549-BX202 


2.4831 


26.2 ± 9.9 


4.2 ± 0.9 


Q1623-BX129 


2.3125 


21.1 ± 2.9 


3.5 ± 0.4 


Q2206-BX151^ 


2.1974 


15.7 ± 0.7 


6.3 ± 0.3 


SSA22a-D17 


3.0851 


25.1 ± 2.5 


2.6 ± 0.3 


SSA22a-C49 


3.1531 


18.0 ± 0.6 


6.4 ± 0.2 


SSA22a-C31 


3.0176 


19.7 ± 1.5 


11.5 ± 0.9 


SSA22a-C32^ 


3.2926 


12.8 ± 35 


5.2 ± 4.7 


Q2343-BX70^ 


2.4086 


16.6 ± 3.5 


4.0 ± 0.6 


Q2343-BX231 


2.4999 


21.5 ± 1.1 


7.2 ± 0.4 


Q2346-BX181 


2.5429 


28.3 ± 2.6 


13.8 ± 1.2 



^ Rest-frame optical nebular emission-line redshifts. Ob- 
jects with z ~ 2.0 — 2.6 have Ho; measurements. Objects 
with z > 3 have [OIII]A5007 measurements. 
^ Observed FWHM and error in units of A. The FWHM 
is a raw value, uncorrected for instrumental broadening. 



^ Emission-line flux and error in units of 10 

cm~^. 
d 



ergs s 



Measured FWHM is less than instrumental resolution. 
^ Measured FWHM is not appreciably higher than the 
instrumental resolution. 

at each wavelength of the true spectrum by a Gaussian 
random number with the standard deviation set by the 
level of the the la error spectrum. Line measurements 
were obtained from the fake spectra in the same manner 
as the actual data. The standard deviation of the distri- 
bution of measurements from the artificial spectra was 
adopted as the error on each centroid, fiux, and FWHM 
value. The NIRSPEC measurements and associated un- 
certainties are given in Table [21 The FWHM values listed 
are not corrected for instrumental broadening. Based on 
Ha or [OHIJA5007 emission line centroid measurements, 
we obtained the systemic redshifts for our sample. These 
redshifts were used to shift both the NIRSPEC and LRIS 
spectra into the systemic rest frame. Based on the Monte 
Carlo uncertainty of the nebular emission line centroid or 
the RMS of the wavelength solution (depending on which 
was greater), we obtain an error of AA ~ 0.9 A, which 
corresponds to a Av ~ 10 — 15 km s~^. Repeat observa- 
tions of two objects from our sample using different slit 
PAs suggest a larger systematic uncertainty on the order 
of 50 km s~^. However, this value should be considered 
an upper limit on the redshift uncertainty given that we 
attempted to match the slit positions between NIRSPEC 
and LRIS observations for a significant fraction of our 
sample. 

4. LYa VELOCITY PROFILES 

Multiple-peaked Lya emission profiles in 2: ^ 2 — 3 
UV-selected galaxies in principle encode a wealth of in- 
formation about the velocity and density fields through 
which they have propagated. It is known that many 
galaxies at this reds hift range are experi encing outfiows 
(jShaplev et "aIll2QQ3l : IVeilleux et al.ll2005[ ). Simple galac- 
tic wind models that incorporate radiative transfer sug- 
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Fig. 3. — Lyo; velocity profiles for the 18 objects in our sample. The zeropoint velocity as measured from the nebular emission (Ho; at 
z r^ 2 or [OIII]A5007 at z ~ 3) is indicated on each image with a dashed line. Considerable variety is evident among the line profile peaks 
and their ratios. From top to bottom and left to right, the panels proceed from Group I to Group V profiles. 



gest that the specific velocity profile of the resultant Lya 
line relative to systemic can be used to determine param- 
eters of the wind such as the expansion velocity, column 
density of windswept material , and outflow temperatures 
(|Verhamme et al.1 I2QQ61 |2QQ8[ ) . With our measured sys- 
temic redshifts and Lyce morphologies, we are now able 
to analyze our systems with this comparison in mind. 
After being shifted to the systemic rest frame based on 
Ha or [OIII] redshifts, the LRIS Lya spectra were trans- 
lated from wavelength to velocity space using the rest 
wavelength of the Lya feature to create the Lya velocity 
profiles. The velocity profiles are shown in Figure [3] and 
the corresponding velocity measurements for each peak 
and trough are shown in Table [3l 

The Lya velocity profiles can be separated into groups 
characterized by the strengths and locations of the emis- 
sion peaks with respect to the velocity- field zeropoint. 
The majority of objects in our sample (11/18), which 
we refer to as Group I, have two peaks that straddle 
the velocity- field zeropoint, with the stronger flux peak 
on the red side. This group includes objects QOIOO- 
MD40, Q0100-BX288, Q0449-BX167, Westphal-DD29, 
Westphal-BX154, Q1549-BX102, Q1549-M22, Q1549- 
BX202, Q2206-BX151, SSA22a-C32, and Q2343-BX70. 
In the case of Q1549-M22, the object with two apparently 
distinct spectral and spatial [OIII] emission components, 
regardless of whether we shift the velocity- field zeropoint 
to Av = 110 km s~^ or Av = — 76 km s~^, as measured 
from the re-extraction of each component individually, 
both Lya peaks still straddle the zeropoint. Therefore 
Q1549-M22 is classified as a Group I object independent 
of the adopted velocity zeropoint. Group II spectra, in- 
cluding objects Q1549-C20 and Q2346-BX181, are de- 
scribed by two peaks that straddle the velocity-field ze- 
ropoint, and have the stronger flux peak on the blue side. 



Group II 



Group III 




Fig. 4. — Two-dimensional Lya emission spectra representing 
each of our 5 groups. For each image the horizontal axis represents 
velocity with respect to the velocity zeropoint, where the dashed 
vertical line represents the systemic velocity (v = 0). The hori- 
zontal axis in each panel spans ±3,000 km s"-*^. The vertical axis 
represents spatial distance along the slit, with a span of ±3.5^^, and 
a vertical scale bar of 1'' is shown in the upper left-hand panel (la- 
beled Group I). Group I (two peaks that straddle the velocity-field 
zeropoint, with the stronger fiux peak on the red side) is repre- 
sented by Westphal-BX154, Group II (two peaks that straddle the 
velocity-field zeropoint, with the stronger fiux peak on the blue 
side) is represented by Q2346-BX181, Group III (both peaks are 
shifted completely redward of the zeropoint, with the stronger fiux 
peak on the blue side) by SSA22a-C49, Group IV (two peaks that 
straddle the velocity-field zeropoint with both peaks having ap- 
proximately the same fiux with the red peak marginally stronger) 
by SSA22a-D17, and Group V (three peaks) by Q1623-BX129. In 
each panel, the locations of the Lya peaks are indicated by arrows. 
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TABLE 3 

Velocity Measurements from Lva Emission and Interstellar Absorption 



Object 


^Hll" 


First Peak^ 


Second Peak^ 


Trough^ 


Absorption^ 


Grism/Grating^ 


Q0100-MD40^ 


2.2506 


-461 


+550 


+288 


-19 


600, 400 


Q0100-BX288 


2.1032 


-579 


+377 


-104 


-317 


400 


Q0449-BX167 


2.3557 


-348 


+586 


-60 


-47 


400 


Westphal-DD29 


3.2401 


-259 


+471 


+202 




300 


Westphal-BX154 


2.5954 


-353 


+404 


-27 


-204 


600 


Q1549-BX102 


2.1932 


-478 


+264 


-155 




600 


Q1549-C20 


3.1174 


-237 


+533 


+210 


+14 


400 


Q1549-M22 


3.1535 


-153 


+414 


+111 




600 


Q1549-BX202 


2.4831 


-542 


+301 


-141 


-143 


600 


Q1623-BX129^ 


2.3125 


-42 


+542 


+229 


-468 


400 


Q2206-BX151 


2.1974 


-417 


+328 


-148 




600 


SSA22a-D17 


3.0851 


-25 


+732 


+385 


-1174 


400 


SSA22a-C49 


3.1531 


+518 


+903 


+735 


-232 


400 


SSA22a-C31 


3.0176 


-15 


+695 


+424 




300 


SSA22a-C32 


3.2926 


-155 


+493 


+99 


-24 


400 


Q2343-BX70 


2.4086 


-476 


+399 


-214 


-33 


400 


Q2343-BX231 


2.4999 


+84 


+567 


+399 




400 


Q2346-BX181 


2.5429 


-214 


+626 


+338 


-49 


400 



^ Rest-frame optical nebular emission-line redshifts. Objects with z ~ 2.0 — 2.6 have Ha measurements. 

Objects with z > 3 have [OIII]A5007 measurements. 

^ Velocity relative to systemic in units of km s"-*^. 

^ All Lya spectra were obtained on the LRIS spectrograph at Keck Observatory. These spectra were 

taken with the 300 line mm~^ grating blazed at 5000 A, or the 400 line mm~^ or 600 line mm~^ grism 

blazed at 3400 A, and 4000 A, respectively, following the LRIS-B upgrade (SteideOraI][200l). 

^ Object Q0100-MD40 had two spectra taken at separate observing runs, with separate gratings. We 

combined both to make the final spectrum. 

^ Object Q1623-BX129 contains a third peak. The corresponding velocity is +1051 km s~^, while the 

second trough is at +846 km s"-*^. 



Objects Q2343-BX231 and SSA22a-C49 comprise Group 
III, where the bluer Lya peak is stronger, but unhke 
the profiles previously described, both peaks are shifted 
completely redward of the zeropoint. Group IV, SSA22a- 
D17 and SSA22a-C31, have two peaks that straddle the 
velocity-field zeropoint with both peaks having approx- 
imately the same flux with the red peak marginally 
stronger. Object Q1623-BX129 accounts for the last 
group. Group V. This object has three peaks. Two of the 
peaks are similar to those of SSA22a-D17 and SSA22a- 
C31, with nearly the same flux and one on either side of 
the zeropoint. The third peak, which is furthest redward 
from the zeropoint, has the smallest flux. Examples of 
two-dimensional Lya spectra from each group are dis- 
played in Figure ID 

A precise velocity-field zeropoint in addition to the ve- 
locity dispersion (see Section 15. ip provides constraints 
on radiative transfer models of the Lya emission from 
our sample of galaxies (see Section [6]). Figure [5] shows a 
graphical summary of peak velocities along with the rel- 
ative fluxes of each peak for every galaxy in our sample. 
Additionally, some objects have detectable interstellar 
absorption lines, which are also labeled on Figure [5l The 
analysis of the absorption lines is described in Section [5^ 
The velocity separation of the peaks in the Lya profiles 
is also an important constraint on possible models. The 
mean velocity peak separation is (Avpeak) = 741 ±39 km 
s~^ (where the error represents the standard deviation of 
the mean), and the median is Avpeak,med = 757 km s~^. 

An examination of the parent LRIS multiple-peaked 
sample indicates the same trends as described above. An- 
alyzing the objects in the parent sample with a S/N lower 
limit of 3 on the strongest flux peak reveals that '^65% 
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Fig. 5. — The location of the velocity for every peak, trough, and 
absorption line with respect to the velocity- field zeropoint. From 
top to bottom, the properties are displayed for Group I, II, III, 
IV, and V objects. Objects are offset vertically for clarity. In each 
line the larger [smaller] black dot represents the velocity peak with 
the larger [smaller] flux. The green square represents the velocity 
of the trough between the peaks. The blue triangle represents 
the absorption line velocity if it was measured for the object. For 
Q1623-BX129, which has three measured peaks, an even smaller 
black dot is used to indicate the third peak with the weakest flux. 
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have their strongest flux peak located on the red side of 
the interpeak trough (compared to 72% from the NIR- 
SPEC sample). Despite a wide range from minimum to 
maximum for the peak-to-peak velocity separation, the 
statistics for the parent LRIS sample are comparable to 
those for our NIRSPEC sample, with a mean value of 
(Avpeak) = 723 ± 18 km s""*^ (where the error repre- 
sents the standard deviation of the mean), and a me- 
dian of Avpeak,med = 721 km s~^. Figure [6] shows the 
distributions of peak-to-peak velocity separation for our 
NIRSPEC sample and the parent LRIS sample. These 
distributions demonstrate that our sample of 18 objects 
with NIRPSEC observations is representative of the par- 
ent LRIS sample as a whole. 

5. PHYSICAL QUANTITIES 

In this section, we discuss several relevant physical 
properties for the NIRSPEC multiple-peaked sample ob- 
jects that can be estimated from our combined dataset. 

5.1. Velocity Dispersion 

As traced by Ha (or [OIIIJA5007) emission, the nebu- 
lar velocity dispersion, a^, reflects the integrated emis- 
sion from star forming regions, which should reveal the 
potential well within the central few kiloparsecs of a 
high-redshift star-forming galaxy. Since the Ha (or 
[OIIIJA5007) velocity dispersion also traces the initial ve- 
locity dispersion of the Lya photons, it is an important 
input parameter for models of Lya radiative transfer and 
will enable a better understanding of the multiple-peaked 
Lya emission systems in our sample. In this paper, aside 
from a brief comparison with previous work (Section l6.2p . 
we do not include a detailed analysis of the relationship 
between intrinsic velocity dispersion and the emergent 
Lya profile. Such analysis is deferred to future work, as 
discussed in Section [71 

The velocity dispersion was calculated as ay = 
(FWHM/2.355) x{c/Xobs), where Xobs is the observed 
wavelength of Ra or [OIII]A5007. The FWHM is deter- 
mined from the subtraction of the instrumental FWHM 
(estimated from the widths of night sky lines) from the 
observed FWHM in quadrature. Uncertainties in the ve- 
locity dispersion were determined using the same method 
described in Erb et al. (2006b). Of our 18 objects, 14 
have measurable velocity dispersion. For 2 additional 
objects we obtained an upper limit in the velocity disper- 
sion because the measured line width was not appreciably 
higher than the instrumental resolution. Finally, 2 ob- 
jects have measured line widths that are smaller than the 
instrumental resolution and are therefore considered un- 
resolved. Figure [7] shows the distribution of the velocity 
dispersion for the sample. The mean value for the sample 
of 14 objects with measurable velocity dispersion is (ay) 
= 90± 9 km s~^ (where the error represents the standard 
deviation of the mean) , while the median for the full sam- 
ple of 18 objects is (Jv,med= 82 kn i s~^. These are si milar 
to the typica l value s reported by lErb et all ()2006bl) and 
iPettini et all (J2001D for z - 2 - 3 UV-selected galaxies 
(cr^ = 108 =b 5 km s-^ for z - 2 and a^ = 84 ± 5 km s"^ 
for z ~ 3). 

5.2. Dust Extinction 



Lya emission can be significantly suppressed due to at- 
tenuation by dust. Commonly, dust extinction in z ~ 2 
star-forming galaxies is estimated fro m rest-fra n ie UV 
colors and an assumption of the Calze tti et al.l ()2000l ) 
starburst attenuation law. Alternatively, the amount of 
dust extinction can be estimated from the observed ratio 
of hydrogen Lyman to Balmer lines, since their intrinsic 
ratios are well described by atomic theory. An assump- 
tion of Case B recombination and T = 10, OOOK deter - 
mines the intrinsic Lya to Ha ratio (Osterbrockl Il989"[ ) . 
Any deviation from this value is then attributed to dust 
extinction. In our sample, only 11 of 18 galaxies have 
measured Ha (the remainder have [OIIIJA5007), making 
the first method more generally applicable for our anal- 
ysis. Both methods were used to calculate E{B — V) for 
objects with measured Ha. 

For the first method, we used the G — 71 colors of our 
sample as determined from Steidel et al. (2003 . 2004) . If 
the Lya line feh in the G-band (at z - 2.48 - 3.28) the 
color was corrected for emission-line contamination using 
Am = 2.5 X log(l + f^), where f^ is the ratio of the 
observed equivalent width of the Lya emission feature to 
the bandwidth of the G filter (1100 A). The correction. 
Am, was then added to the G — IZ color. The G — IZ 
color for objects at z > 2.45 w as also cor rected for IGM 
absorption in the G-band (M adaulll995f ). An intrinsic 
SED model was assum ed as well to es timate the dust 
extinction. We used a Bruzual fc Chari ot (2003) model 
assuming solar metallicity and constant star formation 
rate with a stellar age of 570 Myr for z ~ 2 and 202 
Myr for z ~ 3 objects. Our choice of stell ar ag:e for 
objects at z ~ 2 was based on results from lErb et al] 
(2006b), who found a median age of 570 Myr for star- 
forming galaxies in this redshift range. The stellar age 
of 202 Myr adopted for the z ^3 galaxies was based on 
results from lShaplev et aH (|2005f ). The calculated values 
of E{B - V) are listed in Table El 

For the second method, based on the ratio of Lya and 



Ha emission lines, we adopt the flux ratio 



Lya(int) 



We used the following equation: 



8.7 



Lya^ohs) 



Ha 



(obs) 



= 8.7 X 10"^-"^^'^^"^^'^'^^^""'^^"^ (1) 



In this equation HiLya and HiHa refers to the selective 
extinction at the wavelengths of Lya and Ha, respec - 
tively, as defined in Equation 4 of iCalzetti et al.l (|2000f ). 
E{B — V) can then be calculated based on the observed 
ratio of Lya to Ha flux. All of our calculated E{B — V) 
values from this method exceed those measured from the 
rest-frame UV colors, and the median excess is a factor of 
1.5. This di fference may be explained with reference to 
ICalzetti et al.. (2000.), which states that the stellar con- 
tinuum in local star-forming galaxies should suffer less 
reddening than the ionized gas. Furthermore, Lya may 
suffer additional extinction due to resonant scattering 
and an effectively longer path length though the inter- 
stellar medium. In practice, it is difficult to obtain a 
robust estimate of the Lya to Ha ratio d ue to the effects 
of differential slit loss. lErb et al.l (j2006aD find that NIR- 
SPEC Ha spectra using the .76^^ slit may result in a 
typical factor of ~2 loss, while iSteidel et al.l ()2011l ) find 
that LRIS 1.2^' slit spectra of Lya may result in typical 
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Fig. 6. — Histogram of the distribution of peak-to-peak velocity separation. (Left) The histogram for the 18 objects with NIRSPEC 
measurements. The red shaded bar represents the separation of the 2nd peak (the center peak) from the 3rd peak (the reddest peak) for 
object Q1623-BX129. The mean value for the sample is {Avpeak) = 741 ±39 km s"-*^, with a median of Avpeak,med = 757 km s"-*^. (Right) 
Histogram of the peak-to-peak separation for our parent LRIS sample with a S/N lower limit of 3. The mean value is {Avpeak) = 723zbl8 km 



s ^ , with a median of Avr, 



:721 kms-i. 



losses of a factor of ~5. In detail, the relative loss of Lya 
to Ha depends on the spatial distribution of the emission 
in each of these lines. Overall, the Lya losses appear 
to outweigh those of the Ha losses, which will tend to 
cause an underestimate in the Lya to Ha ratio. Because 
of the systemic uncertainty of the Lya to Ha ratio and 
due to the fact that we can only estimate E{B — V) for 
half of the sample using the (uncertain) intrinsic ratio 
method, we adopt the values calculated from the rest- 
frame UV colors. The average value for our sample when 
employing the first method is {E{B — V))= 0.09, which 
is bluer than the typical value seen a t these redshifts, 
{E{B -V)) = 0.1b (Red dv et all[2008h . This difference 
is not surprising given that our sample is entirely com- 
posed of objects with Lya emission and given the appar- 
ent anti-correlation betwe en Ly a emission EW and dust 
extinction (;Shaplev et aIll20Q3i : iKornei et al.llMlOi ). 

5.3. Star- Formation Rates and Intrinsic Luminosity 

Star-formation rates were calculated from the rest- 
frame UV luminosity. For sources at z ~ 2 we used the 
G-band as a proxy for rest-frame UV luminosity. At z ^ 
3, 7^-band is used as an equivalent probe. The observed 
magnitudes were first corrected for dust extinction using 
rricorr = ^^app — {f^^E{B — V))^ where hz is defined in Sec- 
tion [521 Given the range of redshifts in our z ^ 2 sample 
the G-band effective wavelength of 4780 A corresponds 
to a rest-frame wavelength range of 1330 — 1540 A. For 
the z ^ 3 sample the 7^-band effective wavelength trans- 
lates to a rest-frame effective wavelength of 1590 — 1700 
A. These rest-frame wavelengths were used to determine 
K, values. The corrected magnitude was then used to 
determine the flux, Fjy. Using the luminosity distance, 
we calculated the rest-frame UV luminosity density ac- 




cording to L^ = {AttDl Fu)/{1 + z). 

([1998a,), the relation of SFR(M0 yr-^) = 7.8 x 10" 



29 



100 

CT^ (km s"^) 

Fig. 7. — Histogram of the nebular velocity dispersion, ay, as 
traced by Ho; or [OIII]A5007 measurements. The measurements 
for two objects with FWHM not appreciably larger than the in- 
strumental resolution, Q2206-BX151 and Q2343-BX70, are shaded 
red. The two objects that are considered unresolved, Q0100-MD40 
and SSA22a-C32, are not included in this figure. The mean value 
for the sample of 14 objects with measured velocity dispersions is 
(^av) = 90zb 9 km s"-*^. The median value for the full sample of 18 
objects is CTy^rned = §2 km s"-*^. 



Ljy{eTgs s -^ Hz ^) was used to determine the SFRs as- 
suming a lChabried (|2003[ ) IMF. The SFRs for the sample 
range from 6—110 Mq yr~^. The average value is 32 Mq 
yr~^, which is fairly typical of UV-selected star- forming 
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TABLE 4 
Physical Quantities 



Object 


(Jv ^ 


E{B-V) ^ 


SFR^ 


^Lya 


Q0100-MD40 




0.00 


6 


1.25 


Q0100-BX288 


84.8+^ 


0.02 


11 


2.16 


Q0449-BX167 


50.311^ 


0.20 


51 


10.10 


Westphal-DD29 


102.5+15 


0.00 


73 


14.90 


Westphal-BX154 


^2Atll 


0.26 


110 


21.70 


Q1549-BX102 


64.61^3 


0.10 


13 


2.51 


Q1549-C20 


55.21- 


0.04 


9 


1.83 


Q1549-M22 


153.11^9, 


0.15 


27 


5.23 


Q1549-BX202 


119.7+^^ 


0.07 


14 


2.78 


Q1623-BX129 


89.51^3 


0.00 


9 


1.74 


Q2206-BX151 


< 40.5 


0.01 


8 


1.66 


SSA22a-D17 


12A.^tll 


0.08 


24 


4.84 


SSA22a-C49 


60.5+^ 


0.07 


33 


6.61 


SSA22a-C31 


82.011^ 


0.00 


8 


1.59 


SSA22a-C32 




0.04 


31 


6.19 


Q2343-BX70 


< 78.5 


0.13 


12 


2.35 


Q2343-BX231 


86.3^^ 


0.29 


47 


9.20 


Q2346-BX181 


131. 8+j;^ 


0.18 


97 


19.20 



^ Velocity dispersion measured in units of km s~ . 
^ E{B — V) value derived from G — IZ color. 
*^ Star-formation rate measured in units of M© yr"-*^. 
Intrinsic Lyo; luminosity measured in units of 10 ergs s~ 



galaxies at ^ - 2 - 3 (jErb et aljr2QQ6a[) . The intrinsic 
Lya luminosity, LLyo.^ can be d erived from , the SF R us- 
ing the conversio n described in | Kennicuttl (|1998b[ ) and 
again assuming a iChabrierl ()2QQ3l ) IMF: 



LL2ya(ergs s 



SFR (M0 yr- 
5.1 X 10-43 



(2) 



Based on the average SFR of 32 M© yr"^, the intrin- 
sic Lyof luminosity is (i^Ltya) = 6.3 x 10^^ ergs s"^. The 
intrinsic Lya luminosity is an important quantity for our 
models because it allows us to input an accurate descrip- 
tion of the photon source from our galaxies. SFRs and 
Lhya values are listed in Table ID 

5.4. Interstellar Absorption 

The evidence for feedback in high-redshift galaxies 
comes in several forms, one of which is the nature of 
rest-frame UV interstellar absorption lines. Interstellar 
material is swept up in a galaxy-wide outflow, causing 
a blueshift in the associated absorption lines. In our 
sample over half of the objects have detected interstel- 
lar absorption lines. These include both low-ionization 
(Si IIA1260, Si IIA1526, and C IIA1334) and high- 
ionization (Si IVAA1393,1402 and C IVAA1548,1549) 
lines. 

For all but two objects the absorption lines identified in 
the spectra are low-ionization species. The S/N of these 
lines is not exceptionally strong on an individual basis. In 
principle, all of the low-ionization lines arise in the same 
gas and should share the same Zahs and velocity profiles. 
Based on this idea we continuum-normalized and aver- 
aged the Si IIA1260, Si IIA1526, and C IIA1334 velocity 
profiles to obtain better S/N. Figure [8] shows an exam- 
ple of one of our objects, Westphal-BX154, with aver- 
aged low-ionization lines, as well as object SSA22a-C49, 



which has measurements of both low- and high- ionization 
lines. The velocity offsets in our sample measured from 
the absorption line centroids range from Av^bs = +14 to 
— 1,174 km s~^, relative to the systemic velocity. The 
average value of the interstellar absorption velocity off- 
set for the sample is {Avabs) = — 227 km s~^, which 
is indicative of large scale outflows (iPettini et aDl2QQll : 
iShaplev et aII[2QQ3l : ISteidel et al.ll2QlQ[ r 

The line width of the interstellar absorption features 
was measured as aabs = (FWHM/2.355) x(c/Ao5s), 
where Xobs is the observed wavelength of the interstel- 
lar absorption line. The FWHM was corrected for in- 
strumental resolution {ainst = (FWHMfnst/2.355) = 95, 
155, 190 km s~^ for the 600, 400, 300-line grism, respec- 
tively). The line width, aabs^ is a probe of the range of 
velocities present in the interstellar gas, which is an im- 
portant observational constraint for Lya radiative trans- 
fer models (see Section [6]). The average aabs for the 
12 objects in our sample with absorption line measure- 
ments is 190 ± 63 km s~^ (where the error represents the 
standard deviation of the distribution of measurements). 
This l ine width i s consistent with measurements made 
by Sh aplev et all (|2003 ) of (aabs) = 238 ± 64 km s"^ in 
a composite spe ctrum of 811 z ^ 3 L yman Break Galax- 
ies (LBGs) and iQuider et al.l (|2009f ) of aabs ^ 170-255 
km s~^ in the strong gravitationally lensed galaxy, "the 
Cosmic Horseshoe", at z = 2.38. Our aabs values were 
measured from the combined low-ionization lines for each 
object except for SSA22a-C32, SSA22a-C49 and SSA22a- 
D17, where individual low-ionization lines were measured 
and then averaged. Figure [9] shows the distribution of the 
interstellar absorption line widths for our sample. 

We also investigated if there was an inherent difference 
in the interstellar absorption line profiles for multiple- 
peaked compared to single-peaked Lya emission spectra. 
For this comparison we constructed composite LRIS rest- 
frame UV spectra for our sample of 18 m ultiple-peaked 
objects and for a sample of 29 galaxies from ISteidel et all 
(|201G) with single-peaked Lya emission profiles and pre- 
cisely measured systemic redshifts based on NIRSPEC 
observations of Ha. To create the composite spectrum 
we utilized the IRAF routine, s combine, computing the 
average of the rest-frame spectra at each wavelength with 
a minimum/maximum pixel rejection. The composite 
spectrum was then continuum normalized using a cubic 
spline fitting function. The Lya profile in the compos- 
ite spectrum of our sample is characterized by two clear 
peaks, with a stronger peak on the red side. This pro- 
file is expected since the majority of our profiles resemble 
such a configuration (i.e. Group I). For the single-peaked 
composite spectrum, the measured line widths of the in- 
terstellar absorption lines are systematically larger than 
those in the multiple-peaked composite spectrum, with 
(o'abs) = 140 ± 26 km s~^ for the multiple-peaked com- 
posite and (aabs) = 220 ±19 km s~^ for the single-peaked 
composite. The line widths, (aabs)^ were calculated from 
the averaged Si IIA1260 and C IIA1334 lines, which were 
chosen because of their higher S/N. A Monte Carlo ap- 
proach, similar to the one described in Section [331 was 
used to estimate the error. This difference in interstel- 
lar absorption linewidths may indicate that the range of 
velocities found within the interstellar gas must be nar- 
rower in order to create a multiple-peaked profile instead 
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Fig. 8. — Interstellar absorption velocity profiles for two objects in our sample. Interstellar absorption lines were measured for 12 out of 
the 18 objects in our sample, and the majority have centroids that are blueshifted relative to the velocity zeropoint. (Left) Comparison of 
the Lyo; and low- and high-ionization interstellar absorption velocity profiles for SSA22a-C49. All of the absorption-line velocity centroids 
appear blueward of the velocity- field zeropoint, whereas the Lyo; velocity profile is shifted completely redward of the velocity- field zeropoint. 
(Right) Comparison of the Lyo; and averaged low-ionization interstellar absorption velocity profiles of Westphal-BX154. The low-ionization 
lines. Si IIA1260, Si IIA1526, and C IIA1334, were averaged to obtain better S/N. The absorption line velocity centroid is blueshifted relative 
to the velocity- field zeropoint. The two peaks of the Lyo; velocity profile straddle the velocity- field zeropoint. 
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Fig. 9. — Histogram of interstellar absorption line widths, aabs- 
All aabs values were measured from the combined low-ionization 
lines for each object, except for SSA22a-C32, SSA22a-C49 and 
SSA22a-D17, where individual low-ionization lines were measured 
and then averaged. The mean value of the sample is {(Tabs) = 190 
zb 63 km s~^. The line width is an important parameter for mod- 
eling Lyo; radiative transfer because it reveals the dynamics of the 
absorbing interstellar gas. 



of a single-peaked profile. A comparison of Lya emission 
profiles indicates that the centroid of the Lya emission 
peak in the single-peaked composite is the same as the 
centroid of the redder peak in the multiple-peaked com- 
posite spectrum. Figure [TOl illustrates the differences be- 
tween the total and Group I composite spectra from our 
sample and the s tack of the 29 single-peaked Lya emis- 
sion objects from'Steide l et al.l ([2010). We note here that 
only two out of the 29 single-peaked spectra included in 
the single-peaked composite were taken with the 600-line 
grism, while 600-line spectra comprise 5 out of 18 in the 
multiple-peaked composite. To determine if a possible 
second peak was missed in the single-peaked composite 
due to decreased resolution, we smoothed the multiple- 
peaked composite to the resolution to the single-peaked 
composite. The second peak in the multiple-peaked com- 
posite is still distinctly evident after smoothing, demon- 
strating that the observed multiple-peaked phenomenon 
is not simply a result of better resolution. Also of note 
is the fact that the redshift distribution for our multiple- 
peaked spectra is bimodal (with objects at 2; ^ 2.1 — 2.6 
and z ~ 3.0 — 3.3) and characterized by a mean value 
of (z) = 2.67 ± 0.42, with a median of ^ = 2.54. The 
redshift distribution for the single-peaked spectra has 
{z) = 2.31 ± 0.13, with a median of z = 2.31. To control 
for any possible effects of redshift evolution, we also com- 
pare the subset of z ^ 2 objects with the single-peaked 
spectra, and find the same spectral trends. 

Finally, the interstellar absorption line kinematics al- 
low us to consider the possibility of an alternate sce- 
nario for a multiple-peaked Lya emission profile. The 
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Fig. 10. — Comparison of velocity profiles for the composite of 
the 18 multiple-peaked objects with NIRSPEC measurements, the 
subset of 11 objects that comprise Group I, and 29 single-peaked 
spectra from Steidel et al. (2010). (Top) The Lya velocity pro- 
file for each composite spectrum. The emission centroid from the 
single-peaked composite spectrum matches the centroid of the red 
peak from both of the double-peaked composite spectra. (Bottom) 
The interstellar absorption velocity profile of Si IIA1260 for each 
composite spectrum. The line widths for the single-peaked pro- 
file {{(Tabs) = 220 ± 19 km s"-*^) are systematically larger than for 
the multiple-peaked profile {{(Jabs) = 140 ± 26 km s"-*^). The line 
width for the Group I composite is {(Tabs) — 100 ± 31 km s"-*^. 
The interstellar absorption lines are also inherently deeper for the 
single-peaked compared to the multiple-peaked composite. 



multiple-peaked nature of Lya emission can also be ex- 
plained in terms of the interpeak trough tracing absorp- 
tion from neutral HI with the same kinematics as the 
material giving rise to the interstellar metal absorption 
lines. According to this scenario, the trough and the in- 
terstellar absorption centroids would be aligned in veloc- 
ity space. Figure [5] shows that, with the possible excep- 
tion of Q1549-BX202 and Q0449-BX167, the Lya trough 
velocity centroid does not line up with the interstellar 
absorption centroid for any of the objects in our sample 
with interstellar absorption measurements - ruling out 
this possible hypothesis. 

6. LYa RADIATIVE TRANSFER MODELING 

6.1. Models 

Explaining the emergent profile of Lya emission in 
star-forming galaxies has been the subject of many the- 
oretical and observational studies because of the easy 
observability of the Lya line and the complex processes 
that govern how it emerges from galaxies. Models for 
the propagation and escape of Lya photons have been 
based on both Monte Carlo radiative transfer and sim- 
ple analytic calculations for a var iety of idealized gas 
density and veloci ty distributions (Ha nsen fc Ohll2QQ6l : 
ISteidel et all 1201 Ql : iVerhamme et al.l l2QQ6. 200^ as weh 
as post-processin g of 3D cosmological hydrodynamical 
simulations (.Laursen et al.,, ,2 QQ9a : B arnes et al ,; 201 Ij). 



For analytic models, the gas density and velocity distri- 
butions assumed include static, expanding, and infalling 
shells with both unity and non-unity covering fraction, 
expanding and infalling uniform-density haloes in which 
the gas velocity depends monotonically on galactocen- 
tric distance, open-ended tubes, and spherical distribu- 
tions of discrete clouds with a range of velocity profiles. 
In some of these models, the effec t of dust extinction 
has been treated explicitly (IV erham me et al]|2QQ6l I2QQ8I : 
iHansen fc Ohll2QQ6l : |Laursen et al. 2 QQ9b[ ). 

One of the simple models that has been most closely 
compared with individual LEG spectra is that of an ex- 
pand i ng shell of gaseous material (^Verhamme et al.|[2QQ6l 
I2QQ8I : ISchaerer fc Verhammel[2QQ8[ ). In this picture, en- 
ergy from supernova explosions sweeps up the ISM of 
star-forming galaxies into a geometrically thin but opti- 
cally thick spherical shell. With the advent of efficient 
algorithms to perform resonant-line radiative transfer, it 
has been shown that such a system can be made to repro- 
duce the observe d Lya profiles of high-redshift g alaxies 
(IVerhamme et al.ll2QQ6L boQSD . IVerhamme et al.l (|2QQ6[ ) 
deconstruct the Lya line morphology in order to deter- 
mine how specific features in the Lya profile can be re- 
lated to physical parameters of model systems as well as 
specific photon trajectories. For example, a shell expand- 
ing with uniform velocity ]4xp will produce a primary 
redshifted peak at approximately v = 214xp correspond- 
ing to photons that have undergone a single backscatter- 
ing event en route to the observer. Photons that undergo 
zero backscatterings emerge in an asymmetric double- 



-Ve. 



exp- 



peaked structure centered about a velocity v 
Finally, photons that experience multiple backscattering 
events emerge redward of the primary peak. At fixed 
expansion velocity, the red-side peaks become indistin- 
guishable from one another and the blue-side peak be- 
comes more blueshifted as the Doppler parameter, 6, in- 
creases. Increasing the column density but keeping the 
expansion velocity and thermal velocity fixed results in 
an enhanced inte r-peak separation and broader peaks. 
IVerhamme et al.l (J2006) also consider the effect of dust 
attenuation on the emergent Lya profile, showing how 
it can reduce the overall Lya escape fraction, as well as 
modulate the shape of the profile in velocity space. 

6.2. Comparison Between Models and Observations 

In order to investigate whether such a system can 
reproduce the data presented here, we use the radiative 
transfer code developed in [ Zheng fc Miralda-Escudel 
(2002) and iKollmeier eFall (|2010D to generate 
the observed Lya spectra emerging from expand- 
ing/ infalling shel l models similar to those presented in 
IVerhamme et al.l (|&006), with a range of shell column 
densities, expansion velocities and Doppler parameters. 
Figure [11] shows a subset of such models, for a central 
monochromatic source of Lya photons propagating 
outward through an expanding shell with a column 
density of Nhi = 2 x 10^^ cm~^, expansion velocities 
of Kxp = 100, 300, and 400 km s-\ and Doppler 
6-parameters of 40, 80, 120, and 250 km s~^. It should 
be noted that the Lya radiative transfer calculation 
does not depend on the exact physical dimensions of the 
system for a given column density, expansion velocity, 
and Doppler parameter as long as the shell is kept 
geometrically thin. The models have been smoothed 
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to the typical resolution of our LRIS data, for a fair 
comparison with the observations. Future modeling 
(Kollmeier et al., in prep.) will include a distribution 
of injection frequencies for Lya photons that matches 
the nebular velocity dispersions measured from Ha and 
[OIII] , as well as a more extensive set of gas density and 
velocity distributions for radiative transfer calculations. 



with the typical Group I profile is obtained using the 
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Fig. 11. — Emergent Lya profiles from an expanding shell with a 
central monochromatic source. In all panels, the hydrogen column 
density of the expanding shell is Nhi = 2 x 10^*-' cm~^. The top, 
middle, and bottom panels show the emergent profiles for shells 
with expansion velocities of, respectively, Vexp =100, 300, and 400 
km s~^. In each panel, black, red, green, and blue curves respec- 
tively indicate profiles for Doppler parameters of h =40, 80, 120, 
and 250 km s~^. Each predicted model curve has been smoothed 
to the typical resolution of our LRIS spectra. 



By far the most common Lya profile observed within 
the NIRSPEC multiple-peaked sample is the "Group I" 
type, with two peaks straddling the velocity-field zero- 
point, and a stronger red-side peak. Given the frequency 
of this profile type, it is of particular interest to com- 
pare it with Lya radiative transfer models. In addition 
to the Lya emission properties, we mention here that the 
typical low-ionization interstellar absorption blueshift for 
Group I objects is slightly smaller than the sample av- 
erage, with (^'^o6s) Group I ~ ~^^ ^^ ^~^ (based on 
the 7 out of 11 Group I objects with interstellar absorp- 
tion line measurements). The average low-ionization in- 
terstellar absorption line width for Group I objects is 
(o'abs) ^ 160 km s~^, corresponding to b = 226 km s~^, 
if the line width is interpreted as the thermal velocity 
of a single expanding component of gas (which may in- 
deed not be the correct interpret ation of the inter stel- 
lar absorption profiles, see, e.g., Steidel et al.ll2QlQl and 
Section 16. 3p . As seen in Figure [HI the best match 



model. In this model, which also features two peaks 
straddling the velocity-field zeropoint and a stronger red 
peak, the blue peak appears at Av = —370 km s~^, while 
the red peak appears at Av = +430 km s~^. The model 
matches the average Group I properties both in terms of 
the locations of the Lya peaks relative to zero, as well as 
the interpeak separation, although the predicted contrast 
between the red and blue peak heights is too large. At 
the same time, the predicted absorption profile for such 
an expanding shell (which should be traced well by the 
low-ionization features arising in neutral hydrogen gas) 

would have {(Jabs) ^ b/^/2 ~ 30 km s""*^, i.e., significantly 
lower than what is observed for Group I objects. This 
difference represents a fundamental discrepancy between 
the shell model and the full complement of Group I data 
- even if the average Lya emission profile can be roughly 
reproduced. Figure 



K. 



100 km s- 



other ]4xp = 100 km s" 



/ — ^ /N xw^^ cm ^, 
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smaller column density (i.e., Nhi = 1 x 10 cm~ or 
Nhi = 1 X 10^^ cm~^) and larger Doppler parameters 
{b > 120 km s~^). Overplotted on each panel is the 
continuum-subtracted composite spectrum from Group 
I, which has been normalized to the primary peak height 
of each model. The models with the smaller column den- 
sity also predict emergent profiles with two peaks strad- 
dling the velocity zeropoint and a stronger red-side peak. 
However, in comparison to the Group I composite spec- 
trum, the smallest column density model features a peak 
velocity separation that is too small {Av = 640 km s~^), 
and, in both lower column density cases, the emission 
peaks themselves are too narrow. An important limita- 
tion in this analysis is the oversimplification of comparing 
a model characterized by a specific set of physical param- 
eters with a composite spectrum. While the composite 
spectrum offers a boost in S/N relative to an individual 
spectrum, it also includes the spectra of many individ- 
ual galaxies that likely span a range of parameters (e.g. 
Nni^b.Ve^p). 

The Lya profiles for objects in Group II are also char- 
acterized by two peaks straddling the velocity-field zero- 
point, but with a stronger blue peak. A double-peaked 
profile with a dominant blue pe ak represents a key sig- 
nature of infallin^ g;as ([Zheng^ fc Miralda-Escude 2 0021: 
Diikstra et al.ll2006l : 1Verhamme~et al.ll2006t iBarnes et~al] 



2011D . In terms of the shell models presented in Fig- 
ure [TTJ the predicted Lya profile for an infalling shell 
of material with a specific column density and Doppler 
parameter can be derived by simply fiipping the corre- 
sponding expanding- shell_jnodel about the line Av = 0. 
As sh own inlVerhamme et aTl (J2006) and Diikst ra et al.l 
(2006), infalling spherical haloes of gas will also yield 
double-peaked profiles with a stronger blue peak, assum- 
ing either uniform emissivity or a central point source. 
The lack of significant blueshift in the observed low- 
ionization interstellar absorption lines of Group II ob- 
jects (with Q1549-C20 actually showing a slight red- 
shift) lends additional support to a model of infalling 
gas. In detail, however, there are mismatches between 
the particular infall models mentioned above and the 
observed Group II spectra. The average peak separa- 
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tion for the Group II objects is {Avpeak) ^ 800 km s ^, 
with the blueshifted peak l ocated at Av ^ ~2Q0 km s~^. 
The infahing halo model of IVerhamme et aJ] ()2QQ6l ) pre- 
dicts a peak separation of only (Avpeak) ^ 460 km s~^, 
which is significantly smaller than the observed one (see 
their Figure 5). While an infalling shell model with 
"Onfall = 100 km s~^ and 6 = 40 km s~^ produces roughly 
the correct peak separation and stronger blue peak, the 
locations of the blue and red peaks in the infalling shell 
model relative to zero velocity do not match those of the 
observed Group II spectra. 
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Fig. 12. — Emergent Lya profiles from simple shell models that 
are qualitatively similar to Group I spectra. All panels show 
the predicted profiles for a central monochromatic source passing 
through an expanding shell with Vexp =100 km s"-*^. The top panel 
features a model with Nhi = 2x 10^*-* cm~^and Doppler parameter 
6 = 40 km s"-*^. The middle and bottom panels show the results for 
models with Nhi = 1 x 10^^ cm"^ and Nhi = 1 x 10^'^ cm"^, re- 
spectively. In both of these panels, the predicted profiles are shown 
for Doppler parameters 6 = 120 and 250 km s~^. In all panels, a 
double-peak profile is shown, with two peaks straddling the veloc- 
ity field zeropoint, and a stronger red peak. Overplotted on each 
panel is the Group I continuum-subtracted composite spectrum. 
The Group I spectrum was normalized to match the primary peak 
height of each model. The Nhi = 2 x lO^*-* cm~^ model has the 
closest resemblance to the typical properties of the Group I Lya 
spectra in terms of both precise peak locations and line widths. 
The Nhi = 1 x 10^^ cm~-^ model displays peaks that are both too 
close together in velocity space and too narrow. While the peak 
locations in the Nhi = 1 x 10-*^^ cm~^ model are a fine match to 
the Group I composite spectrum, the predicted line widths are still 
too narrow. 



Like the spectra in Group II, Group III spectra also 
exhibit a stronger blue peak in the Lya profile, but both 
blue and red peaks are redshifted relative to the velocity- 
field zeropoint, and most likely trace outfiowing rather 
than infiowing gas. Only one of the two Group III ob- 
jects (SSA22a-C49) has measurable low-ionization inter- 



stellar absorption. From these features, we measure a 
significant blueshift of Av ~ —230 km s~^, and a line 
width of aabs "^ 230 km s~^. In terms of the simple shell 
models shown in Figure [TTJ the closest match to Group 
III is found for V^xd = 400 km s~ \ and 6 = 40 km s~^. 
As discussed by lYang et al.l ( 120111 ), the blueshifted peak 
in a shell model decreases in strength as Ve^p/b increases. 
This trend is refiected in the large Ve^p/b model discussed 
here, with its lack of significant fiux blue wards of zero, 
and two peaks redshifted relative to zero. Although there 
is qualitative agreement between the model and Group 
III spectra, in so far as both have two peaks redshifted 
relative to zero, and no significant flux blue wards of zero, 
the actual peak separations and locations in the model 
and data do not agree in detail. The model predicts 
^'^peak ^ 600 km s~^, while the Group III objects have 
^Vpeak ^ 400 km s~^. Furthermore, the blue peak in 
the model is shifted by only ~ +100 km s~^ relative to 
zero, while the blue peak in SSA22a-C49 is shifted by 
> +500 km s~^ relative to zero. Finally, the relatively 
small Doppler parameter of 6 = 40 km s~^ would give rise 
to significantly narrower low-ionization interstellar ab- 
sorption line widths than those observed. While Group 
III objects most likely indicate the presence of an outfiow, 
the models described here still present some mismatches 
with the Group III data. 

Distinguishing between Group II and Group III ob- 
jects by establishing the velocity-field zeropoint has po- 
tentially important consequences. Over the last several 
years, the process of galaxy growth through cold gas ac- 
cretion has received much attention in th e theoretical lit- 
eratu re. Both analytic calculations (Birnboim & Dekell 
20031) and numerical simulations (Fardal et al. 2001; 



Kereset al.ll2005L 120091 : iDekel et al.l I2009D suggest that 






high-redshift galaxies primarily grow by smoothly accret- 
ing cold gas from the surrounding IGM. Identifying ob- 
servational signatures of infalling gas is crucial for testing 
the theoretical paradigm of cold accretion. In partic- 
ular, since the predicted signatures of accreting gas in 
low-ionization metal absorption lines ar e so weak, due to 
the lo w metallicity of infalling material (|Fumagalli et al.l 
l2011f ), the Lyof feature may represent the best hope for 
detecting cosmological infall. Beyond the Group II and 
Group III objects in our NIRSPEC sample, - 30% of the 
double-peaked spectra in the LRIS parent sample have a 
stronger blue peak. If the velocity-field zeropoint can be 
established in these objects, to test whether the stronger 
blue peak is also blueshifted relative to zero (i.e., distin- 
guishing between Group II and Group III spectra), their 
Lya spectra may provide evidence for and constraints 
on the nature of gas infall in star-forming galaxies at 
z- 2-3. 

For the remaining three spectra (Groups IV and V), 
we do not find even qualitative matches among the sim- 
ple models presented here. Indeed, for SSA22a-D17 
and SSA22a-C31 (Group IV), there are two peaks with 
roughly equal heights, as seen in a static configuration. 
However, the locations of the peaks are not symmetri- 
cally distributed around the velocity-field zeropoint, as 
would be expected for a static shell, and the interpeak 
trough is redshifted by several hundred km s~^. Further- 
more, SSA22a-D17 has an extremely large low-ionization 
interstellar absorption blueshift measured, with Av ~ 
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— 1170 km s~^, which additionahy rules out the static 
case. Q1623-BX129 has two strong peaks with roughly 
equal heights and similar locations to those of SSA22a- 
D17 and SSA22a-C31, plus a third, redder peak, which 
is significantly weaker. A larger model parameter space 
will be required to match the profile shapes of these three 
objects. 

In summary, while general qualitative matches can be 
found for the Group I, Group II, and Group III Lya 
profile types, we find notable specific discrepancies be- 
tween the models and the totality of the data presented 
here - with particular attention to the interstellar ab- 
sorption profiles. One important issue is that the shell 
models that provide the best matches to, e.g.. Group I 
profiles (Kxp = 100 km s~^, b = 40 km s~^), are not 
able to match simultaneously the broad low-ionization 
interstellar absorption troughs characterized by a me- 
dian of (Tahs.med ^ 190 km s~^. We also highlight the 
importance of our new NIRSPEC measurements for un- 
tangling the nature of these Lya spectra. Measurements 
of the velocity-field zeropoint (i.e. from the Ha or [OIII] 
redshift) are key for distinguishing between Group II (in- 
flow) and Group III (strong outflow) objects. Further- 
more, determining the velocity width of the input Lya 
spectrum (based on the Ha or [OIII] velocity dispersion), 
so that it is not a free parameter of Lya radiative trans- 
fer models, reduces the degeneracies associated with fit- 
ting physical models to Lya profiles. Indeed, when fit- 
ting the double-peaked Lya profile for an object with 
unconstrai ned systemic velo c ity an d input Lya velocity 
dispersion, IVerhamme et al.l (J2008f ) find two expanding- 
shell model solutions that differ by 5 orders of magnitude 
in inferred hydrogen column density (see their Figure 8). 
One of the best-fit models has an inferred intrinsic Lya 
FWHM= 700 km s~^ - a factor of three larger than 
the average for our sample! Furthermore, the placement 
of the velocity-field zeropoint differs b y ^ 500 km s~^ 
for th e two best-fit models in question. IVerhamme et al.l 
(J200o) conclude that neither solution is satisfactory and 
that future observations are needed to constrain their 
models. With systemic redshift and velocity dispersion 
measurements, such degeneracy in best-fit models is no 
longer allowed, potentially leading to much tighter con- 
straints on the underlying physical picture. 

6.3. Additional Considerations and Caveats 

In the preceding discussion, within the framework of 
expanding/infalling shell models, we presented an in- 
terpretation of the observed low-ionization interstellar 
absorption line widths in terms of the random ther- 
mal broadening of gas in the shell. An alternative ex- 
planation of interstellar absorption line widths is pre- 
sented in iSteidel et al.l (|2010l ). according to which ab- 
sorption is produced by a population of discreet clouds 
at a large range of galactocentric radii (i.e., not a thin 
shell). The expansion velocities of these clouds increase 
smoothly and monotonically with increasing radius, and 
their covering fraction is a decreasing function of ra- 
dius. In such a model, in which photons scatter off 
of the TLya = 1 surfaces of discrete clouds, the na- 
ture of the bulk gas kinematics modulates the emer- 
gent Lya profiles, as opposed to the neutral hydrogen 
column density and Doppler 6-parameter of individual 
clouds. This clumpy outflow model accounts for the typ- 



ical interstellar absorption profiles both along the line 
of sight to individual UV-selected star-forming galaxies 
at z ~ 2 — 3, as well as along averaged offset sight- 
lines with impact parameters up to ~ 100 physical kpc, 
and predicts extended Lya emission consistent with the 
observations (S teidel e t al. 2011). At the same time, 
high-resolution 3D cosmological hydrodynamical simu- 
lations that follow the radiative transfer of Lya photons 
through the gas distributions surrounding high-redshift 
galaxies are highlighting the importance of orientation ef- 
fects. The gas column density and velocity distributions 
in these simulated galaxies are anisotropic and irregular, 
characterized by clu mps and filanients of infalling and 
outflowing cold gas (jLaursen et al.l l2009at iBarnes et al.l 
i2011i) . This anisotropy leads to a strong dependence of 
the shape of the emergent Lya profile on viewing an- 
gle. For example, if the line of sight to a galaxy inter- 
cepts a filament of infalling gas, a blue-asymmetric peak 
may be observ ed, even in the pres ence of outflowing gas 
in the system flBarnes et al.ll2Qllf ). A close comparison 
is required between the predictions of such simulations 
and the full ensemble of Lya profiles observed in LBGs, 
as well as incorpo rating detailed Lya radiat i ve tra nsfer 
calculations fro n i IZheng fc Miralda-Escudd ()2002f ) and 
iKoll meier et al.l (120101) into the clumpy outflow model 
of Ste idel etHI (j2010f ). 

Another potentially important effect to consid er is ab- 
sorption by the IGM (Dijkstra et al. 2007; Zhen g et all 
2010). Recent simulations by Laursen et al. (2011) have 
demonstrated that the blue peak in an emergent double- 
peaked Lya profile may be suppressed due to absorption 
by neutral hydrogen external to a galaxy. For the pur- 
poses of such analysis, "IGM" is defined as gas located at 
a radius greater than ~ 1.5 times the dark- matter halo 
virial radius. IGM absorption may be significant even at 
redshifts well below the epoch of reionization. This re- 
sult is relevant to the interpretation of our Group I profile 
types. In such objects, the Lya peak bluewards of the 
systemic velocity is weaker than the red peak. We have 
compared these profiles to simple galactic-scale models 
for outflowing gas, in which the asymmetry arises due 
to the radiative transfer of Lya photons, and have ig- 
nored the possible effects of IGM absorption on the ob- 
served red-to-blue peak ratio. One simple test for the 
importance of IGM absorption consists of examining the 
distribution of red-to-blue peak ratios as a function of 
redshift within the LRIS double-peaked parent sample. 
For this test, we divide the sample into high- {z ^ 3) 
and low-redshift {z ~ 2) subsamples at a threshold red- 
shift of z = 2.7, and calculate the median red-to-blue 
peak ratio for each subsample. If IGM effects are signif- 
icant, we expect to see the red-to-blue peak ratio evolve 
towards lower values at lower redshifts, as the blue peak 
becomes less suppressed by IGM absorption (under the 
assumption that the two samples have the same intrinsic 
distributions of double-peaked Lya profile shapes before 
being subject to IGM absorption). In fact, we find that 
the median red-to-blue peak flux ratios for z ~ 3 and 
z ^ 2 subsamples are 1.3 and 1.5, respectively, and not 
statistically distinguishable. Therefore, IGM absorption 
does not appear to have a significant impact on the ob- 
served multiple-peaked Lya profile shapes in our sample. 

Finally, we must admit the possibility of an entirely 
different explanation for the origin of multiple-peaked 



MULTIPLE-PEAKED LYa EMISSION AT Z - 2 - 3 



17 



Lyof emission profiles, such as galaxy mergers. A 
multiple-peaked profile would naturally result from a 
system in a merger state, with different peaks corre- 
sponding to different m erging components (ICooke et alJ 
120 Id : iRauch et £011201 1). The continuum and emission- 
line morphology of objects in our sample are crucial to 
consider when examining possible merger models. Fig- 
ure [1] shows the (seeing-hmited) 7^-band morphologies 
for objects in our sample. The majority of our tar- 
gets appear to be consistent with single components at 
this resolution. A small number of objects, however, 
are slightly extended (e.g. Q2206-BX151, Q1549-M22) 
or appear to have multiple components (e.g. Q0449- 
BX167, Westphal-BX154). Thre e of th ese objects have 
already been discussed in Section 13.1.11 in terms of their 
complex two-dimensional spectra. If our multiple-peaked 
Lya profiles were simply caused by mergers, there should 
be evidence of multiple-peaked emission in the rest-frame 
optical nebular lines (i.e.. Ha and [OIII]). These features 
more closely trace the systemic velocity and dynamics 
and are typically not as sensitive to bulk outflow mo- 
tions as the Lya emission line. As described in Section 
13.1.11 only Q1549-M22 shows possible signs of nebular 
lines with multiple-peaked behavior in the spectral di- 
rection, ruling out mergers as a possible explanation for 
the majority of our Lya profiles. 

7. SUMMARY AND DISCUSSION 

There is strong observational evidence for the ubiq- 
uity of large-scale outflows in high-redshift star-forming 
galaxies. Gauging the overall impact of these large-scale 
outflows on the evolution of the galaxies sustaining them 
remains an open challenge. Theoretical considerations 
suggest that, at the same cosmic epochs, the inflow of 
cold gas is also a very important process in the evolution 
of star- forming galaxies. Multiple-peaked Lya profiles 
in star-forming UV-selected galaxies at z ~ 2 — 3 of- 
fer a unique probe of both outflows and inflows in the 
early universe. The complex line structure potentially 
provides constraints on many important gas flow param- 
eters. Specifically, the locations of peaks with respect to 
the velocity- field zeropoint, their separations, and flux 
ratios are all additional pieces of information not seen 
in single-peaked Lya profiles. Along with observational 
data such as accurate systemic redshifts, intrinsic nebu- 
lar line widths, and intrinsic ionizing photon fluxes, one 
can distinguish among the possible processes underlying 
these systems. 

We have shown that the phenomenon of multiple- 
peaked profiles appears in a significant fraction (~ 30%) 
of UV-selected star-forming galaxies at 2: ~ 2 — 3 with 
Lya emission. After identifying possible candidates for 
this study we presented a sample of 18 multiple-peaked 
objects with observed Ha or [OIII] nebular emission 
lines, which were used to establish accurate systemic red- 
shift measurements. The average velocity dispersion, ay, 
and the mean star-formation rate in our sample are com- 
parable to values measured from the full population of 
UV-selected star- forming galaxies at z ~ 2 — 3. At the 
same time, given our focus on Lya-emitting galaxies, the 
average E{B — V) in our sample is bluer than what is typ- 
ically measured. Detected interstellar absorption lines 
for objects in the sample show, on average, a blueshift, 
which is suggestive of large-scale outflows driven by su- 



pernovae or winds of massive stars. Additionally, the 
interstellar absorption line widths are measured to be 
{o'abs) ^ 200 km s~^, which is similar to values found in 
other studies of star- forming galaxies at z ~ 2 — 3. A 
comparison of our interstellar absorption line widths to 
a sample of star- forming galaxies with single-peaked Lya 
emission show the single-peaked objects to have consis- 
tently larger interstellar absorption line widths. The dif- 
ference in interstellar absorption line widths may signify 
that the ranges of gas velocities are different for multiple- 
peaked compared to single-peaked Lya emission objects. 

We have qualitatively compared our observed Lya 
profiles with simple models of expanding and infalling 
gaseous shells and halos, given the atte ntion such models 
have recently received in the literat ure ( j Verhamme et al] 
120061 . 120081 : iSchaerer fc Verhammd 120081). The closest 
match between these simple models and our data is found 
for Group I type profiles, which make up 11 out of 18 of 
the NIRSPEC sample objects, and are characterized by 
a red-asymmetric double peak straddling the velocity- 
field zeropoint. Group I type profiles can be produced 
in an expanding shell model, with a high-column den- 
sity (A^ = 2 X 10^^ cm~^), moderate expansion speed 
(Kxp = 100 km s~^), and small Doppler parameter, 
{b = 40 km s~^). At the same time, the observed low- 
ionization interstellar absorption line profiles for Group 
I objects are significantly broader than the features that 
would arise from an expanding shell with 6 = 40 km s~^, 
and signal a problem with the above interpretation for 
the Group I Lya spectra. Certain aspects of Group II 
and Group HI Lya spectra are also reproduced by the 
simple models considered here, though, in detail, the 
matches are not perfect. For the remaining handful of ob- 
jects in our NIRSPEC multiple-peaked sample, the suite 
of models we've considered cannot be used to explain the 
observed profiles, both in terms of Lya peak locations, 
relative strengths, and numbers. 

It is also of interest to compare Lya radiative trans- 
fer models with high-redshift data of significantly higher 
S/N and res o lution than the spectra presented here. 
lOuider et al.l (|2009l ) present the multiple-peaked Lya 
emission profile of the strongly gravitationally lensed 
galaxy, "the Cosmic Horseshoe," for which the magni- 
fication is estimated to be a factor of ~ 24. The Cosmic 
Horseshoe has a very well-defined stellar systemic red- 
shift, from which the Lya spectral profile can be calcu- 
lated iii__vel£cit;Y_s2a£e.^^Q^ an expanding shell model 
from lVerhamme et al] (|2006f ) with Nhi = 7x 10^^ cm"^, 
Vexp = 300 km s~^, and 6 = 40 km s~^, provides a qual- 
itative match to the two redshifted peaks in the Cosmic 
Horseshoe Lya profile, the small Doppler 6-parameter 
in the model is again at odds with the broad low- 
ionization interstellar absorption profile. Furthermore, 
in detail, the locations of the Lya peaks in the predicted 
model spectrum do not line up with those of the obser- 
vations. Both typical and strongly-lensed UV-selected 
star-forming galaxies with well-determined systemic red- 
shifts and, accordingly, Lya profiles calculated robustly 
in velocity space are now providing important obser- 
vational constraints on Lya radiative transfer models. 
These objects highlight the need for considering simulta- 
neously the emergent Lya emission profile, and the pro- 
file of low-ionization interstellar absorption, which must 
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be matched within a unified framework (iSteidel et al.l 

Hoio). 

While we have mainly focused here on the discrep- 
ancies between models and data, and have therefore 
not yet been able to derive the physical parameters of 
the gaseous flows surrounding UV-selected star-forming 
galaxies at high redshift, we view this analysis as an im- 
portant initial step. It will simply not be possible to con- 
strain the gas density and velocity distributions in these 
circumgalactic flows without the simultaneous establish- 
ment of the location of the multiple-peaked Lya pro- 
file and interstellar absorption features in velocity space. 
We have presented such measurements here for a sam- 
ple of objects whose red- and blue-asymmetric multiple- 
peaked Lya profiles at least qualitatively suggest a range 
of processes including both outflows and infall. Addi- 
tional critical observables presented here which will be 
be considered in future modeling are the input velocity 
distribution of Lya photons as traced by the velocity dis- 
persion of either Ha or [OIII] emission, and the degree of 
dust attenuation as traced by rest-frame UV colors and 
Lya/Ra flux ratios. With a realistic model that success- 
fully matches the combination of intrinsic and emergent 



Lyof velocity distributions as well as interstellar absorp- 
tion profiles, we will obtain the gas density and velocity 
distribution as a function of radius. These distributions 
will potentially allow us to constrain the rates at which 
cold gas mass is flowing in and out of galaxies, which 
represents one of the most important goals in the study 
of galaxy formation. 
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